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Abstract 



The classical evolution equations of the Abelian Higgs model are studied at 
temperatures below the Ginsburg temperature of a phase transition which is as- 
sumed to be second order. It is shown that the initial thermal fluctuations provide 
a domain structure which is stable against late time fluctuations. This result lends 
support to the Kibble mechanism for the formation of topological defects. 
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1. Introduction 



Topological defects 1 -* are playing an increasing role in various branches of 
physics. In particular, cosmic strings and global textures give rise to attractive 
scenarios for the origin of structure in the early Universe (for recent reviews see 
e.g., Refs. 2 and 3). It is therefore important to obtain a detailed understand- 
ing of the rate of formation of topological defects in phase transitions from a hot 
symmetric phase to a cold phase with broken symmetry. 

The original mechanism of a defect formation is due to Kibble 1 ). He argued 
that at the phase transition, in any theory which admits topological defects, a 
network of such defects with correlation length (i.e., typical separation) £ will be 
frozen in at the Ginsburg temperature Tq. Here, £ is the correlation length at tQ, 
the cosmic time corresponding to temperature Tq. 

Starting point of Kibble's argument was the assumption that on scales larger 
than £, the orientation of the order parameter in the vacuum manifold is random, 
that the order parameter smoothly interpolates between these random values, and 
that there thus is a finite probability to have nontrivial winding. This probability 
depends on the topology of the vacuum manifold and has been calculated in several 
interesting cases in Ref. 4. 

The Kibble argument has been widely used in cosmology. For example, it has 
been used to generate the initial string configurations 5 ) for cosmic string evolution 
studies, to calculate the abundance of magnetic monopoles 6 ), and to justify the 
occurrence of textures 7 ) in models with nonvanishing 713 (A4), M. being the vacuum 
manifold. Another application of the Kibble mechanism is the prediction of vortices 
in a pressure quench of superfluid helium 8 ). 

The two main assumptions of the Kibble mechanism are that the order param- 
eter takes on random values in M. on scales larger than £, and that it smoothly 
interpolates between its values at different points in space (the "geodesic rule"). 
Recently 9 ), the validity of the geodesic rule has been challenged, in particular 
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for gauge theories. Since the energy density in gauge theories is proportional to 
(-D^0) 2 and not (<9^0) 2 (where is the covariant derivative operator), the or- 
der parameter need not interpolate smoothly in order to minimize the gradient 
energy. Hence, it has become important to investigate the validity of the Kibble 
mechanism more carefully. 

For theories with a global symmetry, the geodesic rule is well justified. In this 
case, the Kibble mechanism has been tested both in numerical simulations 10 -' and 
in the laboratory 11 ' 12 ). 

For theories with a local symmetry the situation is less clear. In models with 
a first order phase transition, a recent analysis 13 ) of the dynamics of the classical 
fields has provided strong support for the Kibble mechanism. The main idea of the 
analysis, however, hinged on the phase transition proceeding via the nucleation and 
subsequent collision of bubbles of the broken symmetry phase. The methods are 
therefore not directly applicable to models with a second order phase transition. 

In this letter we investigate the solution of the classical equations of motion 
for an Abelian Higgs model with a second order phase transition. We consider a 
field configuration set up by thermal fluctuations at the Ginsburg temperature and 
study its stability against thermal fluctuations present at later times. We conclude 
that the initial domain structure is preserved, although naturally the amplitude of 
the order parameter increases. Our results lend support to the hypothesis that the 
Kibble mechanism applies also to gauge theory defects produced in a second order 
transition. 
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2. System and Basic Equations 



As a toy model we consider the Abelian Higgs model with a complex scalar 
field and a U(l) gauge connection A^. Its Lagrangean is 

C = (D^D^ - V(<f>) - \F, V F^ (2.1) 



where 



D„ = - ieA^ (2.2) 



is the covariant derivative with gauge coupling constant e, F^ v is the field strength 
tensor, and V(<p) is the potential for 0. The potential is chosen such that phase 
transition is second order. 

In Lorentz gauge 

= (2.3) 
the equations of motion for and become 

dV 

(d^ - 2ieA^ - e 2 A^)<P + 2 = (2.4) 



and 

d^A v - 2e 2 A v \<j>\ 2 = -ietfdvt . (2.5) 

It is convenient to separate Eq. (2.4) into equations for the amplitude p and phase 
a of 0. Inserting 

= pe ia (2.6) 

into (2.4) we obtain 



8V 

d 2 p - (da - eAfp - e 2 A 2 p + 2-^p = (2.7) 
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and 



d 2 a + 2(d fI a-eA^d^p- = 0. (2.8) 



We wish to study the evolution of the solutions of the coupled system (2.5), 
(2.7) and (2.8) of differential equations during a second order phase transition. We 
assume the following form for the finite temperature effective potential 

Vt(0) = ^A|0| 4 - X -(W - AT 2 )|0| 2 + \\^ . (2.9) 

Here, A is the self-coupling constant of 0, A is a coupling constant determined by 
the graphs which dominate the one loop effective potential, and rj is the value of 
|0| in the vacuum manifold (the scale of symmetry breaking). In the following we 
take A ~ A. T is the temperature. 

The phase transition proceeds by spinodal decomposition 14 ). Due to thermal 
fluctuations in the initial state, the phase a will take on random values on scales 
larger than the correlation length £. The temperature dependence of £ can be 
obtained by equating spatial gradient energy required to set up the domains and 
potential energy gain by having deviate from 0. For temperatures below the 
critical temperature 

T c = A 1//2 A~ 1//2 r7 ~ f] , (2.10) 



the result for f (T) is 1 ) 



/ 2\ - 1 / 2 

^),rVV J • (2.11) 



As long as the temperature is higher than the Ginsburg temperature Tq, these 
domains are unstable to thermal fluctuations. The value of Tq < T c is determined 
by equating thermal energy with the energy of the domains determined above. The 
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result is 



1 



(2.12) 



These qualitative arguments suggest that below Tq, the domain structure is frozen 
in. Since the winding number n 7 around a closed curve 7 is 



the random distribution of phases will induce a finite probability to have a topo- 
logical defect (in our case a cosmic string) in any correlation volume ^(Tq) 3 . 

To reach the final conclusion, we implicitly made use of the "geodesic rule" 
which says that along a curve 7 connecting two domains, the phase a will smoothly 
interpolate between the values of a in the two pluses. This rule can be justified in a 
theory with global symmetry (since it comes from minimizing the spatial gradient 
energy), but not in a gauge theory (since the gauge fields can compensate large 
scalar field gradients). 

In previous work 13 ) , we studied defect formation in a gauge theory with a first 
order phase transition. In this case, the phase transition proceeds by nucleation of 
true vacuum bubbles (|0| ~ rj, a constant) in a surrounding sea of false vacuum 
(0 = 0). Using the equations of motion to study the evolution of <fi when bubbles 
collide, we were able to show that the geodesic rule is dynamically realized. 

In the case of a second order transition we must use rather different methods. 
We will again argue that the geodesic rule is dynamically realized. The basic logic is 
as follows: For T > Tq, thermal fluctuations will dominate the field configurations 
for and A^. We can picture all fields as superpositions of plane waves with 
uncorrelated phases. Fluctuations with wavelength £(T) will dominate. Those 
with £ < £(T) are energetically suppressed while those with £ > £(T) are phase 
space suppressed. 

We will therefore study an initial field configuration with typical wavelength 
£(T) for which the geodesic rule is obviously satisfied. We shall argue - using the 




(2.13) 
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dynamical equations (2.5), (2.7) and (2.8) - that the induced domain structure 
is stable for T < Tq: the change in the phase induced by the dynamics for T < 
Tq is subdominant, and thermal fluctuations at T = Tq do not lead to a large 
perturbation from the evolution of the amplitude p in the absence of fluctuations. 

Note that our analysis is done in a particular gauge. However, the final conclu- 
sion concerning the nonvanishing probability for nonzero winding number is gauge 
independent. 

3. Perturbation Analysis 

In the following we study the initial stages of spinodal decomposition when 
p <C rj. We consider as initial field configuration at T = Tq a plane wave 
Po(x), ao(x), Aq(x) with wave number k (the phases of the three fields are ar- 
bitrary). In the absence of fluctuations this background configuration will evolve 
smoothly in time preserving the initial domain structure and winding number. 

Our main goal is to study the stability of this configuration (and its winding 
number) against thermal fluctuations. The coupling to thermal fluctuations can 
be modelled 15 -* by introducing a second scalar field ip which is assumed to remain 
in thermal equilibrium and which is coupled to via the interaction Lagrangean 

£/ = ^ 2 |0| 2 , (3.1) 

leading to a source term —gi> 2 p on the right hand side of (2.7). 

The effects of thermal fluctuations are analyzed to first order in a perturbative 
expansion about (po, ao, Aq). The total field configurations are 

p(x,t) = p (x,t) +pi(x,t) 

^fe t) = ao(x, t) + ai(x, t) (3.2) 
A(x,t) = Ao(x,t) +A I (x,t) , 

and the equations of motion are expanded to first order in pj, aj, Aj and g. All 
perturbations are induced by the coupling to the thermal bath and hence will be 
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of order g. Our goal is to show that the perturbations remain smaller than the 
background solution provided that k^ 1 corresponds to the Ginsburg length £(Tq). 

The perturbation equations are 



(a? - v 2 ) PJ - 



2 ,„ ,2 dV ddV 
«o 2 -(Vao) 2 -2^ + 2 W p 

[-g^j 2 + do"/ - Va Vaj] Po 



Pi 



(3.3) 



and 



1 1 

(d 2 - V 2 )aj + 2d^a I d flPo - + 2d^a Q d^ Pl - - 2d»a d^ ^ = , (3.4) 

Po PO Po 

where the gauge fields have been set to zero to simplify the equations (they will be 
included later). The perturbations vanish at the initial time t — tQ. 

For k < A 1//2 77 = k cr a the background equation for po has an exponential 
instability 

d 2 p ~ mjpo ~ (A77 2 - k 2 a )po ~ Wpo ( 3 -5) 

with 

Po (x, t) ~ e mi{t - tG \ mi ~ A 1//2 77 . (3.6) 

The solution for the phase «o can be taken to be time independent with spatial 
gradient proportional to k. 

The basic logic of our analysis is as follows: we first analyze the equation of 
motion (3.3) for pj under the assumption that the length scale of fluctuations in 
ai is the same as for those in pj. Using the Green's function method it can be 
shown that 

— < 1 for k p < k crit , (3.7) 
Po 

where k p is the wave number for the perturbation. Next, we consider the equa- 
tion of motion (3.4) for aj and show that modes of aj do not grow in time for 
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wavenumbers larger than k. This shows that the assumption about the length 
scale of fluctuations in aj in the first step of the analysis was self-consistent. In 
a final step we demonstrate that the prescence of gauge fields does not invalidate 
the previous considerations. Now the details: 

Picking out the dominant terms based on the above considerations, Eq. (3.3) 
becomes 

(<9 2 - V 2 )pj - mj Pl = \-gi> 2 - Va Va/]p - ~9^PQ (3.8) 

with the same mj as in the background equation (3.5). This equation can be solved 
using the Green function method, i.e., 

t 

Pl (x, t) = -g J drd 3 yG ret (t-T,x- y)^ 2 (r, y)po(r, y) , (3.9) 


where G re t is the retarded Green function. A new time variable has been chosen 
such that t = is the onset of spinodal decomposition when pj = 0. 

In order to estimate the right hand side of (3.9) we can replace po(r,y) by its 
upper bound 

po(r,y)< A e miT , (3.10) 

where Aq is the initial amplitude. If the phase transition is rapid, we can drop the 
r dependence of the thermal source ip 2 (r,y). Inserting the form of the retarded 
Green function 

1 f e~ ipx 

Gret(x) = -J7T-U / dp- 2 (3.11) 

(27t) 4 J (po + %ey - p 1 + mj 



we obtain 

oo t 



a An f r F p (mi+ipo)T 

Pi^*)~j£ya / fptP'Pip) / dpoe-^ / dr — - . (3.12) 

(27r) 4 J -J J {po + is) 1 ~P l + mi 

-oo o 

The r integral can be done explicitly. At large po, the resulting integral scales as 
Pq 3 . Hence, the integral over po can be written as an integral over a closed contour 
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and solved using the method of residues. It is easy to consider the cases \p\ < mj 
and \p\ > mj separately. In both cases, the integral is bounded by rr/mj. Hence, 
our estimate for (3.12) is 



Pi(t,x) ~ ^-A e mit ^ 2 (x) . (3.13) 

ZH\ 771 j 



The amplitude of ?p 2 (x) is determined by thermal equilibrium. If ^ is a self 
interacting scalar field with self coupling constant of the order 1, then 



ip 2 {x)~T 2 . (3.14) 



Hence, the ratio of pi to po is 



— ~g—- (3.15) 

PO mj 



At the onset of spinodal decomposition 

T = T G <T C = A 1/2 A~ 1/2 r] = A 1/2 ^" 1/2 r] (3.16) 
(see (2.10)), and hence 

^<1. (3.17) 
Po 

We conclude that thermal fluctuations do not dominate the evolution of p for 
T <T C (they do dominate at higher temperatures!) However, the derivation makes 
use of two key assumptions which must be shown to be at least self consistent with 
the above analysis. Firstly, it is assumed that the wave number of fluctuations in 
aj are smaller than k. Otherwise, the second source term in (3.8) (proportional to 
VaoVaj) might dominate. Secondly, it must be shown that the presence of gauge 
fields does not disrupt the smooth evolution of p and a. 
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First, we analyze the evolution of a(t, x). The initial field configuration po(0, x), ao(0, x)| 
is assumed to have nonvanishing winding number along a given circle C in space. 
For winding number 1, the phase cto(0,x) increases smoothly from to 2tt along 
C. The only way the winding can disappear is for a discontinuity in a(x) to de- 
velop. This is only possible if at some point along C the modulus p vanishes. As 
shown above this is not possible as long as ai does not develop small wavelength 
excitations. 

We thus analyze the equation (3.4) for ai under the assumption that pi/po < 1 
in which case the last term on the left hand side of the equation is negligible. We 
approximate the second term as 2ik fi d^aj. In Fourier space, the resulting equation 
is 

(d t 2 + kj - 2k ■ £;/)«/ = 2k ■ kj — (3.18) 

Pi 

where tilde signs denote the Fourier modes labelled by wave vector kj. This inho- 
mogeneous equation can once again be solved by the Green function method. For 
ki ^> k the solution is 

~ (i \ k ' kl Pl( k l) a ^a\ 

a/(fc/,r) -J , (3.19) 

fcj p 

which shown that short wavelength fluctuations are suppressed compared to long 
wavelength inhomogeneities. 

The arguments which lead to (3.17) and (3.19) are self consistent. We could 
also reverse the logic and first study the a equation under the assumption pi/po < 
1, with the result that fluctuations of a on wavelengths smaller than k" 1 are 
suppressed. This result could then be used in the p equation to verify that indeed 
pi/po < 1- 

All the previous arguments neglect the presence of gauge fields and are thus 
only applicable to global symmetries and global defect formation. We now demon- 
strate that the inclusion of gauge fields (working in Coulomb gauge) does not 
invalidate the main conclusions. 
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In the presence of gauge fields, the perturbation equation (3.3) reads 

d 2 P i- 



(da - eA o y + e 2 A 2 - 2 ^ + 2- p 



dV d dV 

|_ 2 ; 

dp 2 dp dp 
[-gi? + e 2 A A T + (da - eA )(d ai - dA T )] p 



Pi 

' (3.20) 



At the beginning of spinodal decomposition (T = Tq), the amplitude of Aq = is 
constrained by thermal equilibrium considerations. Note that the initial configu- 
ration at Tq will typically contain magnetic fields with coherence length £,(Tq) 16 ). 
Hence, it is invalid to set Aq = as initial configuration. However, since 

\D^D^\ < T 4 (3.21) 

and if) ~ T, it follows that A is bounded by 

A<e^T. (3.22) 

Inserting this bound in (3.20) we conclude that the terms involving gauge fields 
are negligible provided 

T < X 1/2 r] = T A . (3.23) 

We suspect that improved estimates would show that this "critical" temperature 
is in fact Tq. 

In the presence of nonvanishing gauge fields, Eq. (3.4) is modified by 

d"ao -> d^a - eA% (3.24) 

and by the presence of an extra source term eA^ However, even at T = Tq the 
new terms are of the same order of magnitude or smaller than the ones present in 
(3.4). Hence, the previous conclusions concerning the absence of short wavelength 
inhomogeneities in aj remain true. 

In summary, we have shown that even in the presence of gauge fields, the initial 
domain structure and winding numbers present at T — Ta are frozen in. 
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4. Discussion 



In this paper we have studied the classical dynamics of spinodal decomposition 
in the Abelian Higgs model coupled to a thermal bath. We considered the classical 
evolution of an initial field configuration at the Ginsburg temperature which was 
assumed to have nontrivial domain structure on a length scale £(7b) (and hence 
a finite probability to have nonvanishing winding number). We argued that the 
initial domain structure is stable against thermal fluctuations. Hence, the winding 
number is preserved. Our results lend support to the Kibble mechanism for defect 
formation - for both global and local defects. 

Above the Ginsburg temperature, thermal fluctuations dominate. Configura- 
tions with nonvanishing winding number are created and destroyed with a fre- 
quency determined by thermal equilibrium. Our results indicate that this domain 
structure created by thermal fluctuations freezes out at Tq, leaving behind a field 
configuration with finite probability for nonvanishing winding number on the scale 
t(T G ). 

Our analysis is unfortunately rather qualitative. It is based on approximate 
equations. We introduced thermal fluctuations by a coupling to an additional 
scalar field, but we did not introduce the dissipation term which must be present 
in order to be in concordance with the fluctuation-dissipation theorem. However, 
such dissipation terms would only slow the growth of fluctuations and would hence 
only strengthen our conclusions. 

The analysis of the equations presented in this paper is only approximate. 
Only the initial stages of spinodal decomposition in which most nonlinearities are 
negligible were analyzed. Even this analysis was done in an approximate way. 

Our classical analysis can be improved by evolving the solutions of the equa- 
tions of motion numerically. A second improvement would be to study the Liouville 
equation associated with our set of classical equations and to derive a result for 
P(tpi, tpi, t), the probability distribution of the fields 15 ). Here, tpi stands for the 
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collection of all fields in the problem. Finally, the analysis should be extended to 
a semiclassical one to control the effects of quantum fluctuations near the onset of 
spinodal decomposition. Work on these projects is in progress. 
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